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We show using computer simulations and mean field theory that a system of particles in two 
dimensions, when confined laterally by a pair of parallel hard walls within a quasi one dimensional 
channel, possesses several anomalous structural and mechanical properties not observed in the bulk. 
Depending on the density p and the distance between the walls Ly, the system shows structural 
characteristics analogous to a weakly modulated liquid, a strongly modulated smectic, a triangular 
solid or a buckled phase. At fixed p, a change in Ly leads to many reentrant discontinuous transitions 
involving changes in the number of layers parallel to the confining walls depending crucially on the 
commensurability of inter-layer spacing with Ly. The sofid shows resistance to elongation but 
not to shear. When strained beyond the elastic limit it fails undergoing plastic deformation but 
surprisingly, as the strain is reversed, the material recovers completely and returns to its original 
undeformed state. We obtain the phase diagram from mean field theory and finite size simulations 
and discuss the effect of fluctuations. 

PACS numbers; 62.25.+g, 64.60.Cn, 64.70.Rh, 61.30.V 



I. INTRODUCTION 



Recent studies on various confined systems [l|, 0, H, 0, 
i, S 0, S S [M [m, E [S M M, M M, M have shown 
the possibiUty of obtaining many different structures and 
phases depending on the range of interactions and com- 
mensurability of a microscopic length scale of the sys- 
tem with the length scale of confinement. Since internal 
structures often play a crucial role in determining local 
dynamical properties like asymmetric diffusion, viscosity 
etc. and phase behavior 0], such studies are of consider- 
able practical interest. Charged particles interacting via 
a screened Coulomb potential when confined in a one di- 
mensional parabolic well showed many zero temperature 
layering transitions At high temperatures this clas- 
sical Wigner crystal melts and the melting temperature 
shows oscillations as a function of density [1]. Such os- 
cillations are characteristic of confined systems, arising 
out of commensurability. In Ref. @, 0] a similar layer- 
ing transition is found in which the number of layers in 
a confined liquid with smectic modulations changes dis- 
cretely as the wall-to-wall separation is increased. The 
force between the walls oscillates as a function of this 
separation depending crucially on the ratio of the sepa- 
ration to the thickness of the smectic layers. Confined 
crystals always align one of the lattice planes along the 
direction of confinement 0, 0] and confining walls gen- 
erate elongational asymmetry in the local density pro- 
file along the walls even for the slightest incommensu- 
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ration [sl]. For long ranged interactions, extreme local- 
ization of wall particles has been observed 0, • Hard 
spheres confined within a two dimensional slit bounded 
by two parallel hard plates show a rich phase behav- 
ior [l|, 13, 0, HO, im, [121 . For separations of one to five 
hard sphere diameters a phase diagram consisting of a 
dazzling array of up to 26 distinct crystal structures has 
been obtained 0, [3| • Similarly confined systems with 
more general interactions also show a similar array of 
phases [l6j. In an earlier experiment on confined steel 
balls in quasi one dimensions (QID), i.e. a two dimen- 
sional system of particles confined in a narrow channel 
by parallel walls (lines), vibrated to simulate the effect of 
temperature, layering transitions, phase coexistence and 
melting was observed A recent study on trapped 

QID solid in contact of its own liquid showed layering 
transitions with increase in trapping potential^lTj. Lay- 
ering fluctuations in such a solid has crucial irnpact on 
the sound and heat transport across the svstem[23l. [2^. 

In this paper, we show that a QID solid strip of length 
Lx confined within parallel, hard, one dimensional walls 
separated by a distance Ly has rather anomalous proper- 
ties. These are quite different from bulk systems in one, 
two or three dimensions as well as from QID solid strips 
with periodic boundary conditions (PBCs). We list be- 
low the main characteristics of a QID confined strip that 
we demonstrate and discuss in this paper: 

1. Re-entrant layer transitions: The nature of the sys- 
tem in QID depends crucially on the density p and 
width of the channel Ly. The number of layers 
of particles depends on the ratio of the interlayer 
spacing (fixed mainly by p) and Ly. With increase 
in channel width Ly at a fixed p we find many re- 
entrant discontinuous transitions involving changes 



2 



in the number of layers parallel to the confining di- 
rection. The phase diagram in the p — Ly plane is 
calculated from Monte Carlo simulations of systems 
with finite size as well as mean field theory (MFT). 
While all the phases show density modulations in 
the direction perpendicular to the wall, we identify 
distinct analouges of the bulk phases viz. modu- 
lated liquid, smectic, triangular solid and buckled 
solid. 

2. Anomalous elastic moduli: A solid characterized 
by a periodic arrangement of particles offers resis- 
tance to elongation as well as shear. The QID con- 
fined solid is shown to have a large Young's modu- 
lus which offers resistance to tensile deformations. 
On the other hand the shear modulus of the sys- 
tem is vanishingly small so that layers of the solid 
parallel to the confinining wall may slide past each 
other without resistance. 

3. Reversible failure: Under externally imposed ten- 
sile strain the deviatoric stress shows an initial lin- 
ear rise up to a limiting value which depends on Ly . 
On further extension the stress rapidly falls to zero 
accompanied by a reduction in the number of solid 
layers parallel to the hard walls by one. However, 
this failure is reversible and the system completely 
recovers the initial structure once the strain is re- 
duced quite unlike a similar solid strip in presence of 
PBCs in both the directions. The critical strain for 
failure by this mechanism decreases with increas- 
ing Ly so that thinner strips are more resistant to 
failure. We show that this reversibility is related 
to the anomalies mentioned above. Namely, the 
confined solid, though possessing local crystalline 
order, retains the ability to flow and regenerate it- 
self. In this manner portions of the QID confined 
solid behave like coalescing liquid droplets. A pre- 
liminary study of this reversible failure mechanism 
was reported in Ref.jl]. 

4. Displacement fluctuation and solid order: Long 
wavelength displacement fluctuations in QID are 
expected to destabilize crystalline order beyond a 
certain length scale. While we always observe this 
predicted growth of fuctuations, in the case of con- 
fined system, the amplitude depends crucially on 
the wall separation. If Ly is incommensurate with 
the interlayer spacing, then local crystalline order 
is destabilized. Otherwise, fiuctuations are kinet- 
ically suppressed in the confined system at high 
densities. Finite size effects also tend to saturate 
the growth of fluctuations. Solid strips of finite 
length therefore exhibit apparent crystalline order 
at high densities both in simulations as well as in 
experiments 0. 

We have used an idealized model solid to illustrate 
these phenomena. Our model solid has particles (disks) 
which interact among themselves only through excluded 



volume or "hard" repulsion. We have reasons to believe, 
however, that for the questions dealt with in this pa- 
per, the detailed nature of the inter particle interactions 
are relatively irrelevant and system behavior is largely 
determined by the nature of confinement and the con- 
straints. Our results may be directly verified in experi- 
ments on sterically stabilized "hard sphere" colloids [2^ 
confined in glass channels and may also be relevant 
for similarly confined atomic systems interacting with 
more complex potentials. Our results should hold, at 
least qualitatively, for systems with fairly steep repulsive 
interactions [H [sl, [13 ■ 

This paper is organized as follows. In the next sec- 
tion, we introduce the model confined solid and discuss 
the geometry and basic definitions of various structural 
and thermodynamic parameters. We then introduce the 
various possible structures with their basic characteris- 
tics in section Hill In section HVl this will be followed by 
the results of computer simulations, in the constant NAT 
(number, area, temperature) ensemble, exploring the de- 
formation and failure properties of this system and the 
relation of the various structures described in section IIIII 
to one another. In section |Vl we provide a finite size 
phase diagram obtained from simulations and compare 
it with an MFT calculation. In section IVII we discuss 
our results with emphasis on the role of long wave length 
fiuctuations in the destruction of crystalline order in low 
dimensions and conclude giving some outlook for future 
work. 



II. THE MODEL AND METHOD 

The bulk system of hard disks where particles i and j, 
in two dimensions, interact with the potential Vij = for 
\rij \ > d and Vij = oo for jr^ | < d, with d the hard disk 
diameter and — rj — the relative position vector 
of the particles, has been extensively [1^ 
studied. Apart from being easily accessible to theo- 
retical treatment [33i experimental systems with nearly 
"hard" interactions |25| are available. The hard disk free 
energy is entirely entropic in origin and the only ther- 
modynamically relevant variable is the number density 
p = N/A or the packing fraction 77 — {TT/A)pd'^. Accu- 
rate computer simulations [Slj of hard disks show that for 
T] > r]f = 0.719 the system exists as a triangular lattice 
which melts below rjm = 0.706. The melting occurs pos- 
sibly through a two step continuous transition from solid 
to liquid via an intervening hexatic phase [Sll [3^. Elas- 
tic constants of bulk hard disks have been calculated in 
simulations [s^. [33 |. The surface free energy of the hard 
disk system in contact with a hard wall has also been 
obtained 35] taking care that the dimensions of the sys- 
tem are compatible with a strain-free triangular lattice. 

The wall- particle interaction potential Kvaii(y) = for 
d/2 < y < Ly — d/2 and — 00 otherwise. Here, evidently, 
Ly is the width of the channel. The length of the channel 
is Lx with Lj; 3> Ly. Periodic boundary conditions are 
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FIG. 1: The confined soUd is shown along with the centered 
rectangular (CR) unit cell. For an unstrained triangular lat- 
tice ax = do and ay — \/iao. Gi, G2 and G3 denote the 
directions of the three reciprocal lattice vectors (RLV). Note 
that the third reciprocal lattice direction G3 is equivalent to 
G2, even in the presence of the walls. 




FIG. 2: (Colour online) Solid: Left panel shows a picture of 
10^ uncorrelated, superimposed configurations of a high den- 
sity ij] = 0.85) solid phase. The wall to wall separation is 
commensurate with a ten layered solid at the given density. 
The colour code is such that red (light) means high local den- 
sity and blue (dark) means low density. The right panel shows 
the corresponding structure factor which shows a pattern typ- 
ical for a two dimensional triangular solid. 



assumed in the x direction(Fig[T]). 

Before we go on to describe the various phases we ob- 
serve in this model system, it is instructive to consider 
how a triangular lattice (the ground state configuration) 
may be accomodated within a pair of straight hard walls. 
For the channel to accommodate ni layers of a homoge- 
neous, triangular lattice with lattice parameter oq of hard 
disks of diameter d, (Figl6|) it is required that, 
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{ni ~ l)ao + d 



(1) 



For a system of constant number of particles and Ly, 
ao is a function of packing fraction 77 alone. We define 
xiViLy) = 1 + 2{Ly — d)/V5ao, so that the above con- 
dition reads x — integer = ni (the commensurate con- 
figuration) and violation of Eq.([T]) implies a rectangular 
strain away from the reference triangular lattice of ni 
layers. The lattice parameters of a centered rectangular 
(CR) unit cell are Ux and Uy (Fig. [1]). In general, for 



a CR lattice with a given Ly 



2{Ly - d)/{ni - 1) 



and Qx = 2/ pay, ignoring vacancies. There are two dis- 
tinct classes of close packed planes in the CR lattice. 
Due to the presence of confinement, even for a triangular 
lattice, the set of planes with reciprocal lattice vector 
(RLV) Gi = perpendicular to the walls are dis- 
tinct from the equivalent set of planes with the RLV's 
G2 = x^cos(l) + y|^sin(|) and G3 = x|f cos(|) - 

^|Esin(|) (FigUI). 

Anticipating some of the discussion in section IIVI we 
point out two different but equivalent "pictures" for 
studying the deformation behavior of confined narrow 
crystalline strips. In the first picture, the stress is re- 
garded as a function of the "external" strain. Using the 
initial triangular solid (packing fraction tjq) as reference, 
the external strain associated with changing at a con- 
stant N and Lj, is e = {L^ — L°)/L°. In this case we 



obtain the stress as an oscillating non-monotonic func- 
tion of e. On the other hand, internally, the solid is, free 
to adjust ni to decrease its energy (strain). Therefore, 
one may, equivalently, calculate strains with respect to a 
reference, distortion-free, triangular lattice at 77. Using 
the definition Sd = e^x - £yy = [ax - ao)/ao - {ay - 
V3ao)/V3ao = cix/ao — ay/VSao and the expressions for 
Qx, tty and oq = 2{Ly — d)/\/3(x — 1) we obtain, 

ni -1 X - 1 



£d = 



X - 1 n; - 1 ■ 



(2) 



where the number of layers n; is the nearest integer to x 
so that Ed has a discontinuity at half -integral values of 
X- For large Ly this discontinuity and Sd itself vanishes 
as 1/Ly for all 77. This "internal" strain Sd is related 
non-linearly to e and may remain small even if e is large. 
The stress is always a monotonic function of Ed- Note 
that a pair of variables ry and Ly (or alternately e and x) 
uniquely fixes the state of the system. 

We have carried out extensive Monte Carlo (MC) sim- 
ulations in the constant NAT ensemble using standard 
Metropolis updates for hard disks (i.e. moves are re- 
jected if they lead to overlaps). The initial condition 
for our simulations is the perfectly commensurate trian- 
gular lattice with a fixed number of layers ri; and at a 
fixed packing fraction 77. After equilibration for typically 
about 10^ Monte Carlo steps (MCS), a strain e is im- 
posed by rescaling Lx- Since the walls are kept fixed, 
this strain reduces 77. The sequence of phases and the 
nature of the transitions among them induced by this 
strain is discussed below. 



III. STRUCTURES AND PHASES 

In QID long wavelength fiuctuationsjs^l are expected 
to destroy all possible order except for those imposed 
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FIG. 3: (Colour online) Buckled phase: A small incommen- 
suration is introduced by increasing the wall to wall sepa- 
ration from the value commensurate to a high density ten 
layered triangular solid at = 0.89. This reduces the pack- 
ing fraction to = 0.85 and produces this phase. The up- 
per panel shows a picture of lO'^ superimposed configurations. 
The colour code for the local densities is the same as Fig[2] 
Note that different portions of triangular solid are displaced 
along y- direction to span the extra space introduced between 
the walls. Lines are drawn to identify this shift in triangular 
regions. The lower panel shows the corresponding structure 
factor where the peak in Gi direction is diminished. Some 
extra weak peaks corresponding to super-lattice reflections 
appear at lower values of the wave-number. 



by explicit (as opposed to spontaneous) breaking of sym- 
metry. The confining potential, in the case of a QID 
system in a hard wall channel, explicitly breaks the 
continuous rotational symmetry down to a Z2 symme- 
try (rotations by angle ±7r). This immediately leads 
to all 2n-atic bond orientational orders, like nematic 
(n = 1), square (n = 2), hexatic (n = 3) etc. which 
remain nonzero throughout the phase diagram. This sit- 
uation is therefore similar to a system in a periodic laser 
potential[38, 39, 40]. Apart from orientational symme- 
try, the confining potential explicitly breaks translational 
symmetry perpendicular to the walls leading to a den- 
sity modulation in that direction. Any fluctuation which 
leads to global changes in the direction of layering or 
in the layer spacing for fixed wall to wall separation is 
strongly suppressed. 

For finite systems (and over finite - but very long - 
times) of confined QID strips one obtains long-lived, 
metastable ordered phases which are observable even in 
experiments [15II . It is these 'phases' that we describe in 
this section. For convenience in nomenclature, we con- 
tinue to use the terms solid, modulated liquid and smec- 
tic to denote these phases, though keeping in mind that 
the distinctions are largely quantitative and not quali- 
tative. For example, a weak solid-like local hexagonal 
modulation is present on top of the smectic layering or- 



der in what we call the smectic phase. The smectic de- 
velops continuously and smoothly from the modulated 
liquid. We also denote the sharp changes in observables 
(eg. strength of diffraction peaks) as "phase transitions" 
though they can never be classified as true phase tran- 
sitions in the sense of equilibrium thermodynamics. We 
show, however, that these phase transitions may be de- 
scribed within MFTs (section |V| although they might 
loose meaning when fluctuations are considered in the 
limit of infinite observation time. 

In our description of the phases we make use of the 
structure factor 

= I] exp(-zG.rjfe)\, 

where — Vj — r; with the position vector of par- 
ticle i. We shall use particularly the values of Sa for 
the reciprocal lattice vectors G = ±Gi(7/) and ±G2(?y). 
Notice that ±G2('7) and ±G3(77) as shown in Fig[T]are 
equivalent directions. A plot of Sg for G = {Gx, Gy) in 
the two dimensional plane gives the diffraction pattern 
observable in scattering experiments. For every phase 
discussed below, the diffraction pattern always shows at 
least two peaks corresponding to ±Gi since they repre- 
sent the density modulation imposed by the walls. The 
relative strengths of these peaks, of course, depend on 
the structure of the phase concerned. 

If the separation between the hard walls is kept com- 
mensurate such that X = n;, an integer, at high den- 
sity we obtain a perfect two dimensional triangular solid 
(Figl2|). The solid shows a diffraction pattern which is 
typical for a two dimensional triangular crystal. We show 
later that appearances can be deceptive, however. This 
triangular "solid" is shown to have zero shear modulus 
which would mean that it can flow without resistance 
along the length of the channel like a liquid. Stretch- 
ing the solid strip lengthwise, on the other hand, costs 
energy and is resisted. The strength of the diffraction 
peaks decreases rapidly with the order of the diffraction. 
In strictly two dimensions this is governed by a non- 
universal exponent dependent on the elastic constants 
In QID this decay should be faster. However, larger 
system sizes and averaging over a large number of con- 
figurations would be required to observe this decay, since 
constraints placed by the hard walls make the system 
slow to equilibrate at high densities. For a general value 
of X the lattice is strained which shows up in the rel- 
ative intensities of the peaks in the diffraction pattern 
corresponding to G2 and G3 . 

A little extra space introduced between the walls start- 
ing from a high density solid phase gives rise to buckling 
instability in y- direction and the system breaks into sev- 
eral triangular sohd regions along the x- direction (Fig [3]). 
Each of these regions slide along close packed planes (cor- 
responding to the reciprocal directions G2 or G3) with 
respect to one another. This produces a buckling wave 
with displacement in y- direction travelling along the 
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FIG. 4: (Colour online) Smectic: A confined triangular solid 
at ?7 = 0.85 is strained in a;-direction to a packing fraction rj — 
0.73. The upper panel shows a picture of 10^ superimposed 
configurations. The colour code for local density is same as 
FiglU The lower left panel shows the density modulation in 
y- direction. The lower right panel shows the corresponding 
structure factor where the only remaining strong peaks are 
in Gi direction identifying a smectic order (solid-like in y- 
direction and liquid- like in other directions). This smectic 
phase, however, possesses hexagonal modulations producing 
weak triangular order leading to faint secondary spots in the 
structure factor. 



length of the solid. In conformity with the quasi two di- 
mensional analog .1.. .21. 22] we call this the buckled solid 
and it interpolates continuously from x = n; to n/ ± 1 
layers. This phase can also occur due to the introduction 
of a compressional strain in a;-direction keeping Ly fixed. 
The diffraction pattern shows a considerable weakening 
of the spots corresponding to planes parallel to the walls 
Sgi together with generation of extra spots at smaller 
wave-number corresponding to the buckled super-lattice. 
The diffraction pattern is therefore almost complemen- 
tary to that of the smectic phase to be discussed below. 
We do not observe the buckled solid at low densities close 
to the freezing transition. Large fluctuations at such den- 
sities lead to creation of bands of the smectic phase within 
a solid eventually causing the solid to melt. 

At low enough densities or high enough in- 
commensuration (x half integral) the elongated density 
profiles in the lattice planes parallel to the walls can over- 
lap to give rise to what we denote as the smectic phase 
(FigU]) in which local density peaks are smeared out in x- 
direction but are clearly separated in y-direction giving 
rise to a layered structure. The diffraction pattern shows 
only two strong spots (5*01) which is typical correspond- 
ing to the symmetry of a smectic phase. We use this fact 



FIG. 5; (Colour online) Modulated liquid: A confined trian- 
gular solid at ?; = 0.85 is strained in i-directon reducing the 
density to 77 = 0.6. The upper panel shows a picture of 10"^ su- 
perimposed configurations. The colour code of local densities 
is the same as Fig[2] The lower left panel shows the density 
modulation in y- direction which is like the smectic phase but 
the modulation dies out at the center. The structure factor 
in the lower right panel shows a ring like pattern which is a 
typical signature of liquid. Superimposed on this are strong 
smectic-like peaks in Gi direction. 

as the defining principle for this phase. Note that unlike 
usual smectics there is no orientational ordering of direc- 
tor vectors of the individual particles since hard disks are 
isotropic. 

At further lower densities the relative displacement 
fluctuations between neighbors diverges and the diffrac- 
tion pattern shows a ring- like feature typical of a liquid 
which appears together with the smectic like peaks in the 
direction perpendicular to the walls. We monitor this us- 
ing the relative Lindemann parameter fs^, given by, 

I = {{u^, - u-,f)lal -t- ((z.^, - vy,f)lal (3) 

where the angular brackets denote averages over conflg- 
urations, i and j are nearest neighbors and m"^ is the a- 
th component of the displacement of particle i from it's 
mean position. This phase is a modulated liquid (Fig[5]). 
The density modulation decays away from the walls and 
for large Ly, the density proflle in the middle of the chan- 
nel becomes uniform. A clear structural difference be- 
tween smectic and modulated liquid is the presence of 
the ring pattern in the structure factor of the latter, a 
characteristic of hquid phase (compare Fig[3]and[5]). We 
must emphasize here again that the distinction between 
a modulated liquid and a smectic is mainly a question 
of degree of the layering modulation. The two structures 
merge continuously into one another as the density is in- 
creased. Also when the modulated liquid co-exists with 
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FIG. 6: Results of NAT ensemble MC simulations of iV = 
rix X Uy — 65 X 10 hard disks confined between two parallel 
hard walls separated by a distance Ly — 9 d. For each 77, the 
system was allowed to run for 10'' MCS and data averaged 
over a further 10® MCS. At 77 = 0.85 we have a strain free 
triangular lattice. Plots show the structure factors SG^,i = 
l(+),2(o) for RLVs Gi(r;), averaged over symmetry related 
directions, as a function of ri. Also plotted in the same graph 
is the Lindemann parameter /(□). The lines in the figure are 
a guide to the eye. 



duces rj. Throughout, 5*02 < •S'gi 7^ 0- This is a conse- 
quence of the hard wall constraint [35| which induces an 
oblate anisotropy in the local density peaks of the solid oS 
from commensuration (nonintegral x). As 77 is decreased 
both 5gi and Sc,^ show a sharp drop at 77 = 77^ where 
X — X* ~ fM — 1/2 [Fig. [7] (inset) ]. For jj < rjc^ we get 
= with Sg^ 7^ signifying a transition from crys- 
talline to smectic like order. The Lindemann parameter 
I remains zero and diverges only below 77 — rjc^ (w rjm) in- 
dicating a finite-size- broadened "melting" of the smectic 
to a modulated liquid phase. 

To understand the mechanical response of the con- 
fined strips, we compute the deviatoric stress built up 
in the system as a function of applied strain. The stress 
tensor of a bulk system interacting via two body cen- 
tral potential has two parts: (i) A kinetic part 
leading to an isotropic pressure and (ii) A virial term 
due to inter-particlc interaction cr^^"*. The free-particle- 
like kinetic component of the stress (3(7^^ — —p5\i, and 
the component due to inter-particle interaction a 



fi^)/S with f^: the v-ih. component of inter- 



particle force, S = LxLy area of the strip. The expres- 
sion for the stress tensor for the bulk system of hard disks 
translates to 1531 



the solid, the layering in the former is particularly strong 
due to proximity effects. 
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IV. MECHANICAL PROPERTIES AND 
FAILURE 



In this section we present the results for the mechan- 
ical behavior and failure of the QID confined solid un- 
der tension. We know that fracture in bulk solids oc- 
curs by the nucleation and growth of cracks [U 0, 
l45j . The interaction of dislocations or zones of plastic 
deformation [4^, with the growing crack tip deter- 
mines the failure mechanism viz. either ductile or brittle 
fracture. Studies of the fracture of single-walled carbon 
nanotubesjl^, show failure driven by bond-breaking 
which produces nano cracks which run along the tube cir- 
cumference leading to brittle fracture. Thin nano-wires 
of Ni are knownji^, [s^l to show ductile failure with ex- 
tensive plastic flow and amorphization. We show that 
the QID confined solid behaves anomalously, quite un- 
like any of the mechanisms mentioned above. It shows 
reversible plastic deformation and failure in the constant 
extension ensemble. The failure occurs by the nucleation 
and growth of smectic regions which occur as distinct 
bands spanning the width of the solid strip. 

We study the effects of strain on the hard disk trian- 
gular solid for a Ly large enough to accommodate only 
a small number of layers n; ~ 9 — 25. The Lindemann 
parameter I diverges at the melting transition [5l|. We 
also compute the quantities 6*0=01 and 6*0=02 ■ 

In Figiniwe show how Sq,^, Sg^ and I vary as a func- 
tion of externally imposed elongational strain that re- 



The presence of walls gives rise to a potential which 
varies only in the y-direction perpendicular to the walls. 
Therefore, strains txx and txy do not lead to any change 
in the wall induced potential. As a consequence the 
conjugate stresses for the confined system a^^ = axx 
and Cxy — cfxy However, a strain eyy does lead to a 
change in potential due to the walls and therefore a new 
term in the expression for conjugate stress appears (ssj. 
^yy = ^yy + with a- = -(E<™> /fy™)/^ where w 
denotes the two confining walls. This expression can be 
easily understood by regarding the two walls as two addi- 
tional particles of infinite massfH^]- Thus, to obtain the 
component of the total stress normal to the walls from 
MC simulations we use 



yy"- g 



^{Ly - yi)5{Ly ~yt~ d/2) 



(5) 



The deviatoric stress, a'- 



ayy, versus strain. 



{eyy = 0) curve for confined hard disks is 
shown in Fig. [7l For f] = r]o (e = 0) the stress due to the 
inter-particle interaction is purely hydrostatic with axx = 
ayy as expected; however, due to the excess pressure from 
the walls the solid is actually always under compression 
along the y direction, thus a^'^ > ayy. At this point 
the system is perfectly commensurate with channel width 
and the local density profiles are circularly symmetric. 
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FIG. 7: (Colour online) A plot of the deviatoric stress a 
versus external strain e obtained from our MC simulations of 
65 X 10 hard disks initially at 77 = 0.85 and Ly = 9d. Data 
is obtained by holding at each strain value for 2 x 10* MCS 
and averaging over a further 3 x 10'' MCS. The entire cycle 
of increasing e(o) and decreasing to zero (□) using typical 
parameters appropriate for an atomic system, corresponds to 
a real frequency of a; « 25K Hz. Results do not essentially 
change for ui = IK Hz — IMHz. The line is a guide to the 
eye. The vertical lines mark the limit of stability of the crystal 
(Cr), the two phase region (Cr+Sm) and the onset of the 
smectic phase (Sm). Inset shows the variation of the critical 
X* with ni, points: simulation data; line: X* ~ 1^1 ~ 1/2- 



Initially, the stress increases linearly, flattening out at 
the onset of plastic behavior at 77 ^ ry^ . At 77^ , with the 
nucleation of smectic bands, cr'^ decreases and eventu- 
ally becomes negative. At rjc^ the smectic phase spans 
the entire system and cr'^ is minimum. On further de- 
crease in 77 towards 77C3, cr*^ approaches from below 
(Fig. [7]) thus forming a loop. Eventually it shows a small 
overshoot, which ultimately goes to zero, from above, as 
the smectic smoothly goes over to a more symmetric liq- 
uid like phase - thus recovering the Pascal's law at low 
enough densities. If the strain is reversed by increasing 
77 back to 770 the entire stress-strain curve is traced back 
with no remnant stress at 7; = 770 showing that the plas- 
tic region is reversible. For the system shown in FigUwe 
obtained 77^ ~ 0.77, 77C2 ~ 0.74 and 77C3 ~ 0.7. As Ly is 

increased, 77^ merges with 7^23 for ni ~ 25. If instead, 
and Ly are both rescaled to keep x = "^i fixed or PBCs 
are imposed in both x and y directions, the transitions 
in the various quantities occur approximately simultane- 
ously as expected in the bulk system. Varying Ux in the 
range 50 — 5000 produces no qualitative change in most 
of the results. 

Reversible plasticity in the confined narrow strip is in 
stark contrast with the mechanical response of a simi- 
lar strip in the absence of confinement. In order to show 
this, we study a similar narrow strip of 65 x 10 hard disks 
but now we use PBCs in both the directions. At pack- 



FIG. 8: (Colour online) A plot of the deviatoric stress cr 
versus external strain e obtained from our MC simulations of 
65 X 10 hard disks, in presence of PBCs in both the directions, 
initially a triangular lattice at 77 = 0.85. Data is obtained by 
holding at each strain value for 2 x 10'* MCS and averaging 
over a further 3 x lO** MCS. The entire cycle of increasing e(o) 
and decreasing to zero (□) using typical parameters appropri- 
ate for an atomic system, corresponds to a real frequency of 
ijj ~ lOOK Hz. The line through the points is guide to the eye. 

ing fraction 77 = 0.85 the initial geometry (i^ — 65ax, 
Ly — lOuy with Qy — ^/Sux/'^) contaius a perfect tri- 
angular lattice. We impose strain in a fashion similar 
to that described in FiglT] The resulting stress-strain 
curve is shown in Fig[5] With increase in strain e the 
system first shows a linear (Hookean) response in the de- 
viatoric stress cr = axx — f^yy , fiattening out at the onset 
of plastic deformation below e ^ 0.1. Near e = 0.1 with 
a solid-liquid transition (solid order parameter drops to 
zero with divergence of Lindeman parameter at the same 
strain value) the deviatoric stress a decreases sharply to 
zero in the liquid phase obeying the Pascal's law. Unlike 
in the confined strip, with further increase in strain a 
does not become negative and fluctuates around zero in 
absence of wall induced density modulations. With de- 
crease in strain, starting from the liquid phase, the stress 
cr does not trace back its path. Instead first a large nega- 
tive stress is built up in the system as we decrease strain 
up to e '--^ 0.04. With further decrease in e, the stress 
cr starts to increase and at e = the system is under a 
huge residual stress l3acP = 25. The configuration of the 
strip at this point shows a solid with lattice planes ro- 
tated with respect to the initial stress-free lattice. This 
solid contains defects. Note that in the presence of hard 
confining walls, global rotation of the lattice planes cost 
large amounts of energy and would be completely sup- 
pressed. The generation of defects is also difficult in a 
confined system unless they have certain special charac- 
teristics which we describe later. 

For a confined QID strip in the density range 77^3 < 
77 < 77ci we observe that the smectic order appears within 
narrow bands (Fig. [5]) . Inside these bands the number of 
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FIG. 9: (Colour online) Plot of as a function of the 

channel length x / d a.t ri — 0.76 after time t = (a)5 x 10^ and 
(b)2 X 10*^ MCS for = 10^ x 10 with Ly = 9d. Note that x = 
10 in the solid and x = 9 in the smectic regions. Arrows show 
the coalescence of two smectic bands as a function of time. 
Calculated diffraction patterns for the (c) 10-layered solid and 
(d) 9-layered smectic regions, (e) Close up view of a crystal- 
smectic interface from 10"^ superimposed configurations at 77 = 
0.77. The colour code of the local density is the same as 
Fig[2] Note the misfit dislocation in the inter-facial region. 
A Burger's circuit is shown to identify the Burger's vector 
corresponding to this dislocation as 6 = yay/2 + xax/2. 



layers is less by one and the system in this range of rj is 
in a mixed phase. A plot (FiglS] (a) and (b)) of x(a^ji), 
where we treat x as a space and time (MCS) depen- 
dent "order parameter" (configuration averaged number 
of layers over a window in x and i), shows bands in which 
X is less by one compared to the crystalline regions. After 
nucleation, narrow bands coalesce to form wider bands 
over very large time scales. The total size of such bands 
grow as 77 is decreased. Calculated diffraction patterns 
(Fig. [9] (c) and (d)) show that, locally, within a smectic 
band Sg^ 3> Sg^ contrast to the solid region where 
Sgi ~ Sg2 7^ 0- 

Noting that x — X* — ni — 1/2 when the solid fails 
(Fig. [7] inset), it follows from Eq. [21 the critical strain 



(4n,-5)/(2n,-3)(2n, 



This is sup- 



ported by our simulation data over the range 9 < rt; < 14. 
This shows that thinner strips (smaller n/) are more re- 
sistant to failure. At these strains the solid generates 
bands consisting of regions with one less particle layer. 
Within these bands adjacent local density peaks of the 
particles overlap in the x direction producing a smec- 
tic. Within a simple density functional argument [s^ 



it can be shown that the spread of local density profile 
along X-axis ax ^ 1/ \J Sg^ ^'Hd that along y-direction is 
ay - 1/^25gi + Sg2 [61|. In the fimit Sg^ ^ (melt- 
ing of solid) ax diverges though ay remains finite as Sgi 
remains positive definite. Thus the resulting structure 
has a smectic symmetry. 

A superposition of many particle positions near a solid- 
smectic interface [see Fig. IHe)] shows that: (1) The 
width of the interface is large, spanning about 10 — 15 
particle spacings. (2) The interface between n; layered 
crystal and n; — 1 layered smectic contains a dislocation 
with Burger's vector in the y- direction which makes up 
for the difference in the number of layers. Each band 
of width s is therefore held in place by a dislocation- 
anti-dislocation pair (Fig . In analogy with classical 
nucleation theory [sl, l56l | the free energy Ft of a single 
band can be written as 



-dFs + 



ao 



(6) 



where if'^ is an elastic constant, b = yay/2 + xax/2 is 
the Burger's vector, 5F the free energy difference be- 
tween the crystal and the smectic per unit length and Ec 
the core energy for a dislocation pair. Bands form when 
dislocation pairs separated by s > ■^b'^ / 5F arise due 
to random fluctuations. To produce a dislocation pair 
a large energy barrier of core energy Ec has to be over- 
come. Though even for very small strains Sd the elastic 
free energy becomes unstable the random fluctuations 
can not overcome this large energy barrier within finite 
time scales thereby suppressing the production of n; — 1 
layered smectic bands up to the point of e*^. In principle, 
if one could wait for truly infinite times the fluctuations 
can produce such dislocation pairs for any non-zero 
though the probability for such productions exp(— /?i?c) 
[(3 = l/ksT, inverse temperature] are indeed very low. 
Using a procedure similar to that used in Ref. [sl, 0, [HI] , 
we have monitored the dislocation probability as a func- 
tion of 77 (Fig fTU)) . For confined hard disks, there are 
essentially three kinds of dislocations with Burger's vec- 
tors parallel to the three possible bonds in a triangu- 
lar solid. Only dislocations with Burger's vectors having 
components perpendicular to the walls, cause a change in 
ni and are therefore relevant. The dislocation formation 
probability is obtained by performing a simulation where 
the local connectivity of bonds in the solid is not allowed 
to change while an audit is performed of the number of 
moves which tend to do so. Since each possible distortion 
of the unit cell (see FigHU]- inset) can be obtained by two 
specific sets of dislocations, the dislocation probabilities 
may be easily obtained from the measured probability of 
bond breaking moves. Not surprisingly, the probability 
of obtaining dislocation pairs with the relevant Burger's 
vector increases dramatically as 77 ^ rjc^ (Fig fTU)) and ar- 
tificially removing configurations with such dislocations 
suppresses the transition completely. Band coalescence 
occurs by diffusion aided dislocation "climb" which at 
high density implies slow kinetics. The amount of im- 
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FIG. 10: Dislocation probabilities of a 65 x 10 system are 
plotted as a function of strain starting from a triangular solid 
at 77 = 0.85 and Ly = 9d. The corresponding bond-breaking 
moves are depicted in the inset panel. Arrows show the direc- 
tions of the bond-breaking moves. -I- show dislocation prob- 
abilities for those Burger's vectors which have components 
perpendicular to the walls (topmost figure in the inset). On 
the other hand probabilities corresponding to the other two 
types of moves (denoted by *) remain zero. 



posed strain fixes the total amount of solid and smectic 
regions within the strip. The coarsening of the smec- 
tic bands within the solid background in presence of this 
conservation, leads to an even slower dynamics than non- 
conserving diffusion. Therefore the size of smectic band 
L{t) scales as L ~ i^/'^[53l. Throughout the two-phase 
region, the crystal is in compression and the smectic in 
tension along the y direction so that a is completely de- 
termined by the amount of the co-existing phases. Also 
the walls ensure that orientation relationships between 
the two phases are preserved throughout. As the amount 
of solid or smectic in the system is entirely governed by 
the strain value e the amount of stress a*" is completely 
determined by the value of strain at all times regardless 
of deformation history. This explains the reversible [58j 
plastic deformation in Fig. \7\ 



V. THE MEAN FIELD PHASE DIAGRAM 

In this section we obtain the phase diagram from 
Monte Carlo simulations of 65 x n; hard disks, n; be- 
ing the number of layers contained within the channel. 
We compare this with a MFT calculation. At a given 
Ly we start with the largest number of possible layers of 
particles and squeeze up to the limit of overlap. Then 
we run the system over 5 x 10'' MCS, collect data over 
further 5 x 10^ MCS and increase in steps such that 
it reduces the density by Si] = 0.001 in each step. At a 
given Ly and n; the value of L^ that supports a perfect 



VSir ni{ni - 1) 

8 Ly{Ly ~ 1) 



(7) 



In this channel geometry, at phase coexistence, the stress 
component Uxx of the coexisting phases becomes equal. 
The other component Gyy is always balanced by the pres- 
sure from the confining walls. Thus we focus on the quan- 
tity Vx = —f^xx as Lx is varied: 



Vx 



1 dF 
Ly dLx 



= P 



> dfN dfN _ Px 

dp ' dp 



(8) 



The free energy per unit volume /y — {N/V){F/N) = 
pIn- One can obtain the free energy at any given density 
p by integrating the above differential equation starting 
from a known free energy at some density po, 



fv{p) = PfN{Po) + P 



Px, 

—dp 

PO P 



(9) 



We discuss below how f%{po) for solid and modu- 
lated liquid phases are obtained. The free energy of 
the solid phase may be obtained within a simple ana- 
lytical treatment viz. fixed neighbor free volume theory 
(FNFVT), which we outline below. More detailed treat- 
ment is available in Ref.[27]. The free volume Vf{ri,x) 
may be obtained using straight forward, though rather 
tedious, geometrical considerations and the free energy 
f%{ri, x) — ~plii''^/('7j x)- The free volume available to a 
particle is computed by considering a single disk moving 
in a fixed cage formed by its nearest neighbors which are 
assumed to be held in their average positions. The free 
volume available to this central particle (in units of cP) 
is given entirely by the lattice separations b — ao(l-l-ea;a:) 
and h — •\/3ao(l + ^yy)/'^ where oq is lattice parameter 
of a triangular lattice at any given packing fraction i] 
and e^x = («; - ^ 1), fyy = ^/^xx- As stated in 

Secini X is obtained from channel width Ly and packing 
fraction rj. Vf is the area available to the central test 
particle. Note that the effect of the confining geome- 
try is incorporated in the lattice separations b, h. The 
FNFVT free energy has minima at all x = n;. For half 
integral values of x the homogeneous crystal is locally un- 
stable. Although FNFVT fails also at these points, this 
is irrelevant as the system goes through a phase transi- 
tion before such strains are realized. In the high density 
triangular solid phase, we know f%{pQ = Pa), exactly, 
from the fixed neighbor free volume theory (FNFVT) . It 
is interesting to note that, apart from at p/\, this FN- 
FVT can estimate the solid free energy quite accurately 
at small strains exx < 4% around pA- To obtain the 
solid branch of free energy fy{p) from simulations, we 
integrate the Px — p curve starting from the FNFVT es- 
timate at p = Po = p A • 

For the confined fluid phase we use a phenomenologi- 
cal free energy [63| of hard disk fluid added with a sim- 
ple surface tension term due to the scaled particle the- 
ory [6^. We outline this in the following. For the bulk 
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(unconfined) liquid phase one can obtain f%{p = pi) up 
to a fairly large volume fraction rji ~ 0.81 where the 
liquid reaches its metastable limit in bulk, using the phe- 
nomenological expression [6^. 

(277, - 1) ln(l - "^^^m) - ln(l - ^) 

PInKPi) - TTf, N 

2(1 - ?7c) 

+ lnpi-1 (10) 

where r^c is the close packed volume and depends on Ly. 
We have checked from the simulations that this equa- 
tion of state shows good agreement with the simulated 
Px — V diagram in a confined channel up to a density 
of ?7 = 0.65. Above this density, falls below and py 
goes above this equation of state estimate. Integrating 
from r]Q — rji = 0.6 using the above mentioned form of 
f%{pi) we find the bulk part of the liquid branch of free 
energy fy [p) . In order to incorporate the effect of con- 
finement approximately, we add a scaled particle theory 
(SPT) [g^l surface energy contribution to the bulk free 
energy. Within SPT, the interfacial tension (interfacial 
energy per unit length) of a hard disk fluid in contact 
with a hard wall can be calculated as 7 = 7 + psptd/2 
[13 where 7 = {pd/2(3)[l/{l - 77) - 1/(1 - r])^]a.nd the 
SPT equation of state is Ppspt/p = 1/(1 - ?7^) [sl. Thus 
the SPT surface energy per unit volume St = 1 / Ly = 
pd/2(3{l — r])Ly is added with the bulk contribution fy{p) 
to account for the confinement. 

At thermodynamic equilibrium, the chemical potential 
p = —{dfv/dp)v.T and pressure for the coexisting phases 
must be same. Thus from a plot of free energy vs. density 
obtaining common tangent to solid and liquid branches 
of free energies means finding coexisting densities of solid 
and liquid having the same chemical potential. Once the 
two branches of free energy are obtained in the above 
mentioned way, we fit each branch using polynomials up 
to quadratic order in density. The common tangent must 
share the same slope and intercept. These two conditions 
uniquely determine the coexisting densities. The onset of 
buckling instability can be found from the slope of the 
Px — V curve at high densities - a negative slope being 
the signature of the transition. The estimates of the co- 
existing densities and onset of the buckling instabilities in 
our simulations are thus obtained and shown using filled 
symbols in the phase diagram fPig fTT]) . 

To obtain the phase diagram from MPT (lines in 
Fig fTTj) we use the above mentioned phenomenological 
free energy of hard disk fluid [63| together with the SPT 
surface tension term. Por the solid we use the PNPVT 
free energy for all densities. We take care to minimize 
the free energy of the solid with respect to choices of ni 
first. Then we find out the minimum of the solid and 
liquid free energies at all densities. The coexisting den- 
sities are obtained from common tangent constructions 
for this minimum free energy. Our MPT phase diagram 
so obtained is shown along with the simulation results in 
the 77 — Ly plane fPig fTTj) . 

In FigUT] points are obtained from simulations and the 
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FIG. 11: (Colour online) Phase diagram. The continuous lines 
denote the phase boundaries obtained from mean field theory. 
The white, vertically shaded and horizontally shaded regions 
indicate single phase, two-phase coexistence and forbidden re- 
gions, respectively. The single phase regions of low densities 
are liquid and high densities are solid phases. On the green 
(dashed) lines ni — 9, 10, 11 the system is commensurate with 
the channel width such that exactly integral number of layers 
are contained. Points denoted by filled o, filled A, and filled 
□ are obtained from simulations. Regions with densities less 
than the filled o are modulated liquid and above filled A are 
solid. In the intervening region, we have a solid-smectic co- 
existence. The filled □ mark the onset of buckling instability. 
Smectic order increases smoothly and continuously as density 
increases. The liquid within the coexistence region is highly 
modulated. 



continuous curves from theory. The regions with den- 
sities lower than the points denoted by filled circles are 
single phase modulated liquid. Whereas the regions with 
densities larger than the points denoted by filled trian- 
gles indicate the single phase solid. The modulations in 
the liquid increase with increasing density and at high 
enough p the structure factor shows only two peaks typ- 
ical of the smectic phase. This transition is a smooth 
crossover. All the regions in between the filled circles 
and filled triangles correspond to solid-smectic coexis- 
tence. The filled squares denote the onset of buckling at 
high densities. The high density regions shaded by hor- 
izontal lines are physically inaccessible due to the hard 
disk overlap. The MPT prediction for the solid-liquid 
coexistence is shown by the regions between the continu- 
ous (red) lines (shaded vertically). The unshaded white 
regions bounded by red lines and labeled by the number 
of layers denote the solid phase. All other unshaded re- 
gions denote liquid phase. Por a constant channel width, 
a discontinuous transition from liquid to solid via phase 
coexistence occurs with increasing density. However, the 
MPT predicts that the solid remelts at further higher 
densities. Notice that the simulated points for the on- 
set of buckling lie very close to this remelting curve. 
Since the MPT, as an input, has only two possibilities 



of solid and fluid phases, the high density remelting line 
as obtained from the MFT may be interpreted as the 
onset of instability (buckling) in the high density solid. 
The MFT prediction of the solid-fluid coexistence region 
shows qualitative agreement with simulation results for 
solid-smectic coexistence. The area of phase digram cor- 
responding to the solid phase as obtained from simulation 
is smaller than that is predicted by the MFT calculation. 
This may be due to the inability of MFT to capture the 
effect of fluctuations. From the simulated phase diagram 
it is clear that if one fixes the density and channel width 
in a solid phase and then increases the channel width 
keeping the density fixed, one finds a series of phase tran- 
sitions from a solid to a smectic to another solid having a 
larger number of layers. These re-entrant transitions are 
due to the oscillatory commensurability ratio x that one 
encounters on changing the channel width. This is purely 
a finite size effect due to the confinement. It is important 
to note that for a bulk hard disk system, solid phase is 
stable at 77 > 0.719, whereas, even for a commensurate 
confined strip of hard disks, the solid phase is stable only 
above a density 77 = 0.75. With increase in incommen- 
suration the density of the onset of solid phase increases 
further. This means that confinement of the QID strip 
by planar walls has an overall destabilizing effect on the 
solid phase. 

The phase diagram calculated in this section is a MFT 
phase diagram where the effects of long wavelength (and 
long time scale) fluctuations arc ignored. For long enough 
QID strips fluctuations in the displacement variable 
should increase linearly destroying all possible sponta- 
neous order and leading to a single disordered fluid phase 
which is layered in response to the externally imposed 
hard wall potential. However, even in that limit the lay- 
ering transitions from lower to higher number of layers 
as a function of increasing channel width, survives^. 

Our simple theory therefore predicts a discontinuous 
solid-fluid transition via a coexistence region with change 
in density or channel width. However, details like the 
density modulation, effects of asymmetry in density pro- 
file, vanishing displacement modes at the walls and most 
importantly nucleation and dynamics of misfit disloca- 
tions crucial to generate the smectic band mediated fail- 
ure observed in simulations are beyond the scope of this 
theory. Also, the effect of kinetic constraints which stabi- 
lize the solid phase well inside the two phase coexistence 
region is not captured in this approach. We believe, nev- 
ertheless, that this equilibrium calculation may be used 
as a basis for computations of reaction rates for address- 
ing dynamical questions in the future. 

Before we end this section, it is useful to compare our 
results with similar calculations in higher dimensions viz 
hard spheres confined within hard parallel plates forming 
a quasi two dimensional filmp]. [3. Y2^. Commensuration 
effects produce non-monotonic convergence of the phase 
boundaries to that of the bulk system. The appearence 
of novel phases e.g. the buckled solid not observed in 
the bulk is also a feature showed by all these systems; 
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FIG. 12; (Colour online) Shear stress vs. shear strain at 
77 = 0.85. A system of 40 x 10 hard disks simulated with 
periodic boundary conditions and Ly commensurate with ten 
layered triangular solid gives a shear modulus /i = 398±4. On 
the other hand when the same system of 40 x 10 hard disks is 
confined within a commensurate channel, that fits 10-layers of 
lattice planes, the shear modulus drops drastically to ^ = 0. 

on the other hand these quasi two dimensional systems 
show two kinds of solids - triangular and square - and no 
'smectic like' phase. The effect of fluctuations also should 
be much less in higher dimensions and one expects the 
observed phases to have less anomalous properties. 



VI. DISCUSSION AND CONCLUSION 

One of the key definitions of a solid states that a solid, 
as opposed to a liquid, is a substance which can retain its 
shape due to its nonzero shear modulus. Going by this 
definition, a QID solid strip even of finite length confined 
within planar, structureless walls is not a solid despite its 
rather striking triangular crystalline order. Indeed, the 
shear modulus of the confined solid at 77 = 0.85 is zero, 
though the corresponding system with PBC show finite 
nonzero shear modulus (Sec Fig fT^ . This is a curious 
result and is generally true for all values of 4 < 71; < 25 
and 77 investigated by us. Confinement induces strong 
layering which effectively decouples the solid layers close 
to the walls allowing them to slide past each other rather 
like reentrant laser induced melting [H, [H, |43|. This im- 
mediately shows that the only thermodynamically stable 
phase in confined QID channel is a modulated liquid, the 
density modulation coming from an explicit breaking of 
the translational symmetry by the confinement. 

To understand the nature and amount of fluctu- 
ations in the confined QID solid we calculate the 
auto-correlation of the displacement along the channel, 
B{x) — {{u^{x) — u^(0))^) for a layer of particles near 
a boundary. The nature of the equilibrium displacement 
correlations ultimately determines the decay of the peak 
amplitudes of the structure factor and the value of the 
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FIG. 13: (Colour online) Flucuation < (u^(a;)-u^(0))2 > /d^ 
on a lattice plane (line) for 5000 x 10 commensurate triangu- 
lar solid {ri = 0.85) with PBC in both directions is compared 
with that in presence of hard channel confinement. Averages 
are done over 100 configurations, o denotes the data obtained 
in presence of PBC in both directions and □ denotes the data 
obtained for confined system. In both cases the initial log- 
arithmic growth (tracked by dashed and solid straight lines 
respectively) crosses over to a faster linear increase before 
saturating due to finite size effects. Confinement apparently 
reduces the amount of fluctuations without afecting the gen- 
eral behavior. 



equilibrium elastic moduli [55|. It is known from har- 
monic theories, in one dimension B{x) ~ x and in two di- 
mensions B{x) ~ ln(x). In the QID system it is expected 
that for small distances, displacement fluctuations will 
grow logarithmically with distance which will crossover 
to a linear growth at large distances with a crossover 
length Xc ^ Ly In Ly [s^, [sj] . We calculate this correla- 
tion for a system of 5000 x 10 particles averaged over 
10, 30, 50, 100 configurations each separated by lO'^ 
MCS and equilibrated over 10^ MCS at = 0.85. We 
compare the results obtained for a strip with confine- 
ment and a strip with PBC's in both directions, taking 
care that in each case the channel width is commensurate 
with the inter-lattice plane separation. With an increase 
in the number of configurations over which the averages 
are done, fluctuations reduce and converge to a small 
number. This happens for both the cases we study. It is 
interesting to notice a logarithmic to linear cross-over of 
the fluctuations near x = lOd for both the cases. Since 
the harmonic theoryjs^l ignores effects due to bound- 
aries, one can conclude that this general behavior is quite 
robust. At large distances, displacement correlations arc 
expected to saturate due to finite size effects. The mag- 
nitude of fluctuation at saturation for the confined sys- 
tem 0.0035) is noticably lower than that in presence 
of PBC (- 0.0045) (FigdHD. Thus the introduction of 
commensurate confinement can have a stabilizing effect 
of reducing the overall amount of fluctuations present in 
a QID strip with PBC. 

We have separately calculated the elastic modulus of a 



40 X 10 conflned hard disk solid at = 0.85 under elon- 
gation (Young's modulus) in the x and the y directions. 
For an isotropic triangular crystal in 2D these should be 
identical. We, however, obtain the values 1361 and 1503 
(in units of ksT /cP and within an error of 3%) respec- 
tively for the two cases. The Young modulus for the 
longitudinal elongation is smaller than that in the direc- 
tion transverse to the confinement and both these values 
are larger than the Young modulus of the system under 
PBC {Y = 1350). This corroborates the fact that the 
non-hydrodynamic component of the stress Uxx — <yyy is 
non-zero even for vanishingly small strains as shown in 
FigEl Therefore even if we choose to regard this QID 
solid as a liquid, it is quite anomalous since it trivially 
violates Pascal's law which states that the stress tensor in 
a liquid is always proportional to the identity. This is be- 
cause of the explicit breaking of translational symmetry 
by the confinement. Lastly, commensurability seems to 
affect strongly the nature and magnitude of the displace- 
ment fluctuations which increase dramatically as the sys- 
tem is made incommensurate. 

As Lx ~ Ly ^ oo all such anomalous behavior is 
expected to be localized in a region close to the walls 
such that in the central region a bulk 2d solid is recov- 
ered. This crossover is continuous, though oscillatory 
with commensurabilities playing an important role, and 
extremely slow (~ 1/Ly). It is therefore difficult to ob- 
serve in simulations. 

What impact, if any, do our results have for realis- 
tic systems? With recent advances in the field of nano 
science and technology [1^ [6^ new possibilities of build- 
ing machines made of small assemblage of atoms and 
molecules are emerging. This requires a knowledge of 
the structures and mechanical properties of systems up 
to atomic scales. A priori, at such small scales, there 
is no reason for macroscopic continuum elasticity theory, 
to be valid [67t. Our results corroborate such expecta- 
tions. We have shown that small systems often show 
entirely new behavior if constraints are imposed lead- 
ing to confinement in one or more directions. We have 
also speculated on applications of reversible failure as ac- 
curate strain transducers or strain induced electrical or 
thermal switching devices [2^, [l^l • We believe that many 
of our results may have applications in tribologyflGj in 
low dimensional systems. The effect of corrugations of 
the walls on the properties of the confined system is an 
interesting direction of future study. The destruction of 
long ranged solid like order should be observable in nano 
wires and tubes and may lead to fluctuations in transport 
quantities [68j . 
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